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CHAPTER  1 


INTRODUCTION 

In  this  thesis  a  comparison  of  simulated  and  predicted 
performance  of  a  complex  crosscorrelatlon  receiver  against  a  range 
spread  scatterer  Is  discussed.  Performance  Is  measured  by  the  ratio 
of  the  receiver  output  signal  energy  to  receiver  output  for  noise 
only.  This  ratio  can  be  predicted  assuming  analytic  expressions  for  a 
singly  spread  target  scattering  function  and  the  transmit  and  receiver 
waveforms.  The  slgnal-to-nolse  ratio,  at  the  output  of  the  receiver 
SNR,  Is  examined  as  a  function  of  receiver  Integration  time,  for  a 
fixed  Input  SNR.  The  prediction  equation  was  developed  by  Ricker,  and 
Is  derived  In  terms  of  the  range  spread  function  and  the  receiver 
crossambiguity  function  [1].  The  simulation  Is  a  digital 
Implementation  of  the  assumed  target  and  signal  models. 

Two  variations  of  the  specific  scattering  model  are  considered. 
The  first  target  consists  of  both  a  low  amplitude  continuous 
scattering  component,  and  two  large  amplitude  point  reflectors.  In 
this  model,  the  low  level  continuous  scatterer  Is  represented  as  a 
sequence  of  closely  spaced  point  reflectors,  whose  reflection 
characteristics  are  random  samples  from  uncorrelated  Gaussian 
processes  [2].  A  second  scattering  model  adds  varying  degrees  of 
correlation  to  the  continuous  scatterer  by  filtering  the  uncorrelated 
sequence.  The  degree  of  correlation  Is  inversely  proportional  to  the 
bandwidth  of  the  filter  used,  while  the  actual  shape  of  the  filter 


response  Is  relatively  unimportant.  SNR  comparisons  between  the 
simulated  correlated  scattering  values  and  the  predicted  levels, 
assuming  uncorrelated  scattering,  for  several  filter  bandwldths  will 
be  discussed. 

Chapter  2  presents  the  analytic  model  for  the  range  spread 
scattering  function  and  receiver.  The  amplitude  of  the  scattered 
signal  Is  the  convolution  of  the  transmit  pulse  with  the  Impulse 
response  of  the  random  target.  For  the  thesis.  It  Is  assumed  the 
medium  may  be  modeled  as  a  single  deterministic  propagation  path  with 
any  randomness  In  the  received  signal  due  solely  to  the  scattering 
process.  The  specific  target  model  Is  characterized  by  a  range  spread 
scatterer  [3].  The  median  reflection  characteristics  of  the  scatterer 
are  specified  by  the  range  scattering  function,  S(t),  and  Is  a 
function  of  delay  only.  The  scattering  function  to  be  examined 
consists  of  two  Independent  components,  a  low  level  continuous 
scatterer  and  a  set  of  distinct,  high  amplitude  point  reflectors. 

The  response  of  the  receiver  to  this  signal  Is  measured  by  the 
crossambiguity  function,  defined  as  the  magnitude  squared  of  the 
crosscorrelatlon  between  transmit  and  receiver  processing  signals.  To 
simplify  the  form  of  the  ambiguity  function,  both  signals  are  assumed 
to  be  cosine  pulses  of  frequency  «c>  The  transmit  pulse  length  Is 
Tjf  and  the  processing  length  T,  where  It  Is  assumed 
T  >  T$.  Further,  It  Is  assumed  target  scattering  process  does  not 
spread  the  signal  In  frequency,  and  for  convenience,  the  scatterer 


velocity  Is  assumed  to  be  zero.  At  zero  Doppler  and  delay,  the 
cross-ambiguity  function  reduces  to  simply  the  ratio  of  the  transmit 
pulse  length  to  the  receiver  Integration  time;  l.e.,  T  /T. 

Clearly,  If  the  scattering  function  were  an  Ideal  point 
reflector,  the  received  signal  would  be  simply  a  time  delayed  version 
of  the  transmit  pulse,  and  the  optimum  receiver  Integration  time 
corresponds  to  T  -  T$  [4].  When  the  received  signal  Is  spread  In 
time  by  the  scattering  process,  this  receiver  Integration  time  Is 
suboptimum,  or  "mismatched."  The  performance,  as  measured  by  the 
slgnal-to-nolse  at  the  receiver  output,  will  Increase  for  some  longer 
Integration  time,  T,  that  best  matches  the  received  signal  duration. 

The  simulation  of  the  scattered  signal  and  receiver  are 
discussed  In  Chapter  3.  Both  models  are  derived  from  VanTrees  [5]. 
The  scattering  function  Is  obtained  by  combining  the  models  for  a 
range  spread  scatterer  and  a  slowly  fluctuating  point  target.  The 
continuous  scattering  component  Is  modeled  as  many  closely,  and 
equally,  spaced  point  reflectors.  The  amplitudes  of  these  reflectors 
are  uncorrelated,  complex  Gaussian  numbers  of  equal  variance.  The 
point  reflectors  are  modeled  by  two  uncorrelated,  complex  6auss1an 
numbers,  with  large  variances  compared  to  the  continuous  scatterers. 
For  convenience,  the  point  reflectors  are  placed  along  the  scatterer 
so  that  the  signal  reflected  from  these  points  are  separated  by  at 
least  a  transmit  pulse  length,  l.e.,  the  reflections  from  the  point 
scatterers  are  resolved  by  the  transmit  pulse. 


The  receiver  model  used  Is  shown  In  Figure  1.  The  reception,  of 
either  the  reflected  signal,  or  noise.  Is  multiplied  by  the  complex 
conjugate  of  the  receiver  processing  signal,  g(t).  Integrated  and 
squared.  It  Is  assumed  the  signal  and  noise  are  uncorrelated, 
allowing  the  slgnal-to-nolse  ratio  to  be  formed  by  passing  the  signal 
and  noise  through  the  receiver  separately.  By  averaging  many 
simulation  outputs  for  signal  or  noise,  a  comparison  statistic  for  the 
ambiguity  function  analysis  may  be  generated. 

A  comparison  of  the  simulated  SNR,  and  that  predicted  by  the 
ambiguity  function  analysis  Is  presented  In  Chapter  4.  The  three 
combinations  of  the  scattering  model  are  1)  continuous  scatterer,  2) 
composite  scatterer,  with  both  the  continuous  scatterer  and  point 
reflectors,  and  3)  point  reflectors  only.  All  comparisons  are 
presented  as  a  function  of  the  ratio  of  the  processing  time,  T,  to 
transmit  pulse  length,  T  .  The  agreement  between  simulation  and 
prediction  Is  quite  good,  with  an  average  difference  of  roughly  0.5 
dB,  Independent  of  the  specific  type  of  scatterer,  or  the  value  of 

T/V 

Chapter  5  examines  the  effects  of  adding  correlation  to  the 
continuous  scattering  component.  Correlation  Is  added  by  filtering 
the  uncorrelated  continuous  scattering  sequence,  with  the  filter 
bandwidth  Inversely  proportional  to  the  degree  of  correlation  In  the 
filtered  sequence.  Comparing  the  predicted  SNR  with  simulated  values. 


the  quality  of  comparison  decreases  with  decreasing  filter  bandwidth, 
or  Increasing  correlation  In  the  continuous  scattering  component.  The 
correlated  scatterer  Increases  the  received  signal  energy  since  the 
scattered  signal  now  contains  a  coherent  component  as  well  as  the 
random  scatter.  For  very  narrow  bandwldths,  the  coherent  return  Is 
large,  and  the  simulated  SNR  Is  greater  than  the  prediction  for  all 
values  of  T/T$. 

The  effects  of  correlation  on  the  continuous  scatterer  In  the 
composite  target  Is  also  compared  with  predicted  SNR,  assuming 
uncorrelated  scattering.  For  a  small  degree  of  correlation,  the 
comparison  between  simulated  and  prediction  Is  affected  only  for  data 
prior  to  the  first  point  scatterer  return.  As  with  the  continuous 
scatterer  only  comparisons,  the  simulated  values  are  higher  than  the 
prediction,  but  this  effect  Is  overshadowed  by  the  point  scatterer 
returns  for  relatively  wide  filter  bandwldths.  For  the  very  narrow 
bandwidth,  the  correlated,  or  coherent,  scatterer  energy  obscures  even 
the  point  scatterer  returns,  with  differences  between  simulation  and 
prediction  as  great  as  17  dB.  As  with  the  correlated  continuous 
scatterer  only,  the  added  energy  due  to  the  correlated  scatterers  Is 
Independent  of  T/Ts- 

A  summary  of  the  major  sections  of  this  thesis  Is  presented  In 
Chapter  6.  A  number  of  areas  for  additional  research  are  also 


discussed. 


CHAPTER  2 


MODEL  ANALYSIS  USING  SCATTERING  AND  AMBIGUITY  FUNCTIONS 

2.1  Introduction 


This  chapter  develops  the  analytic  models  for  a  range  spread 
scatterer  and  complex  correlation  receiver.  By  specifying  the  form  of 
the  scattering  function,  transmitted  pulse  and  receiver  processing 
signal,  the  expected  performance  of  the  receiver  may  be  evaluated 
analytically  In  terms  of  the  scattering  function  and  receiver 
ambiguity  function.  Performance  Is  measured  In  terms  of  the  receiver 
output  SNR,  for  a  given  transmit  pulse  length  T  ,  and  receiver 
processing  time,  T. 

2.2  Complex  Envelopes  and  Narrowband  Signals 


The  waveforms  considered  In  this  thesis  are  restricted  to  the  real 
signals  that  satisfy  the  relation 


s(t) 


Re 


s(t)e 


jwct 


(2.1) 


where  ?(t)  Is  the  complex  envelope  of  the  signal,  s(t),  and  w.  is 
the  signal  carrier  frequency  In  radians/second.  These  signals  have 
Fourier  transforms  centered  about  t  «c,  and  are  referred  to  as 
narrowband  [6]. 


A  property  of  narrowband  signals  Is  that  the  complex  envelope  Is 


a  slowly  varying  function  compared  to  the  carrier  frequency.  It  Is 
convenient  to  use  the  complex  envelope  for  linear  filter  operations, 
and  then  take  the  real  part  of  the  output  to  recover  the  real  output 
signal  [7].  The  narrowband  signals  of  concern  In  this  thesis  are  the 
transmit  signal,  f(t),  and  the  receiver  processing  signal,  g(t). 

Further  It  Is  assumed  that  both  signals  are  pulsed  cosines  of  the  form; 


f(t)  -  g(t)  -  Re 


Ae 


0  <  t  <  T1 
elsewhere 


(2.2) 


where  A  represents  the  constant  signal  amplitude.  These  signals  are 
energy  normalized  by  assuming 


dt  -  1 


Combining  Equations  2.2  and  2.3  yields 


f(t) 


2.3  The  Range  Scattering  Function 


0  <  t  <  T1 


(2.3) 


(2.4) 


2.3.1  Definition.  The  received  signal  is  assumed  to  be  the 
convolution  of  the  transmit  waveform  f(t),  with  one  realization  of  the 
time  varying  Impulse  response.  While  In  general,  the  scattering 


process  may  disperse  the  transmit  signal  In  both  time  and  frequency, 
the  scatterer  Is  considered  to  be  stationary  In  this  thesis.  The 
reflected  signal  is  delayed  and  temporally  dispersed  without  any 
Doppler,  or  frequency  spread.  This  type  of  scattering  process  Is 
referred  to  as  range-spread  [8]. 


A  physical  model  of  the  range-spread  scatterer  might  be  a  line 
segment  whose  surface  Is  random  and  rough  compared  to  the  wavelength 
of  the  transmit  signal.  The  signal  reflected  from  any  Increment,  &L, 


along  the  scatterer  Is 

s(t)  f(t-r)  b(  at)  (2.5) 

where  ET  Is  the  energy  In  the  transmit  pulse,  and  ?(t-  t)  Is  the 
complex  envelope  of  the  transmit  signal,  delayed  byr.  b(Ar)  Is  a 
sample  function  from  a  random  reflection  process  and  ir«  2&L/c,  and  c 
Is  the  speed  of  sound. 


The  projected  extent  of  the  scattering  function  Is  assumed  to  be 
much  greater  than  the  transmit  pulse  length  and  Independent  of  aspect 
angle.  The  reflected  signal  Is  then  of  greater  duration  than  the 
transmit  pulse,  or  spread  In  range  by  the  target.  The  reflected 
signal  Is  the  superposition  from  many  locations  along  the  target.  The 
signal  at  some  time  t  Is  written. 


s(t)  £  ^(t-T^)  (t^  atj 


(2.6) 


where  at^  represents  the  time  delay  from  the  AL^  Increment  along 
the  target.  It  Is  assumed  that  the  number  of  scatterers  along  any 
Increment  Is  large  so  that  'B(t)  Is  a  sample  function  from  a  complex 
Gaussian  process  [9]. 


The  statistics  of  the  received  signal  are  completely  specified  by 
the  mean  value  and  the  covariance  function  of  the  random  range  spread 
scatterer  Tj(t) .  The  covariance  function  of  the  reflection  process, 

K  (T.T-j)  Is  defined  as  [10] 


-  E  b(r)b  (t,) 


(2.7) 


The  model  Is  simplified  further  by  assuming  the  reflection  process  Is 
zero  mean  [11],  and  spatially  uncorrelated.  Equation  2.7  may  then  be 


written  as, 


KCr.^)  -  .  E  b(r) 


(2-8) 


The  range  scattering  function  may  now  be  defined  as  [12], 


Sr(t)  «  E  b(r) 


(2.9) 


The  correlation  function  of  the  reflected  signal  in  the  absence 


of  noise  Is, 


Mt,t,)  »  E  s(t)s  (t.) 


(2.10) 
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for  a  scattered  signal  of  duration  Ty  where  Tj  >  T$.  The  total 
energy  In  the  signal  reflected  from  a  zero  Doppler  scatterer  Is  [13] 


E«  *  jf  srs<T-T1>dT  ■  E 

From  Equation  2.13,  note  that  the  scattering  function  at  any  time 
delay, t  ,  represents  the  average  energy  reflected  from  a  portion  of 
the  target  at  range,  r  *  c  f/2. 
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SR(T)dT 


(2.1 


2.3.2  Composite  Range  Scattering  Function.  The  particular 
scattering  function  of  Interest  consists  of  two  Independent,  random 
components.  The  first  Is  a  continuous  scatter,  characterized  by  a 
scattering  density,  p The  scattering  function  for  the  continuous 
component  Is  a  constant,  or 


S  (t) 
CSy  1 


(2.1 


where  p E  Is  In  units  1/seconds.  The  second  scattering  process  Is  a 
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where  represents  the  energy  scattered  from  the  I**1  point 
scatterer.  The  total  target  scattering  function  Is  then. 


SR(t)  "  PE  +  £  Ei*(T-Ti> 
1-1 


(2.15) 


(2.16) 


Figure  2  Is  an  example  plot  of  the  composite  target  scattering 
function  for  comparatively  large  point  scatterer  amplitudes  and  a  low 
continuous  scatter  amplitude.  Convolution  of  the  scattering  function 
and  a  10  millisecond  transmit  pulse  yields  the  reflected  signal  shown 
at  the  bottom  of  the  figure. 


2.4  Receiver  Ambiguity  Function 


2.4.1  Receiver  Model .  The  receiver  used  to  process  the  signal 
reflected  from  the  composite  scatterer  Is  shown  In  Figure  1,  and  Is 
defined  by  VanTrees  as  the  complex  correlation  receiver  [14].  For 
detection,  the  reflected  signal  Is  multiplied  by  the  conjugate  of  the 
receiver  processing  waveform,  g(t).  Integrated  for  a  time  T,  magnitude 
squared,  and  the  output  compared  with  a  threshold. 


FIGURE  2.  Magnitude  of  the  composite  scattering 

function  (top),  and  reflected  signal  (bottom) 
Note  the  difference  between  the  time  axes. 


The  performance  of  the  receiver  is  calculated  as  a  function  of 


the  Integration  time  to  pulse  length  ratio,  T/Ts  >  1,  and  Is 
measured  by  the  ratio  of  the  mean  signal  energy  to  mean  noise  energy 
at  the  receiver  output.  The  output  SNR  Is  then  defined  as 


SHR  ■  tUttl/tUl 
E{n(t)n  (t)| 


(2.17) 


Noise  Is  assumed  to  be  a  white,  band  limited,  6auss1an  noise  of 
power  spectral  density  Ng.  The  output  noise  energy  Is 


E»  ■  "o  /T|9(t,f  dt  *  *0  • 

which  Is  numerically  equal  to  the  the  spectral  density  by  Equation  2.3. 


(2.18) 


2.4.2  Receiver  Ambiguity  Function.  A  measure  of  the  receiver 
response  to  a  returned  signal  Is  the  crossambiguity  function  defined  as 


X(T.«)  2  - 
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f  f(t-f)g*(t+f)eJwt  dt 

mo 

(2.19) 


Note  that  the  term 


*0 

X(t.«)  -  / 7(t-|)g*(t+f)  ejwt  dt 


(2.20) 


Is  the  crosscorrelation  of  the  transmit  and  processing  signals,  and 
Indicates  the  "match"  between  the  reflected  signal  and  transmit  signal 


For  the  assumed  transmit  signal  and  receiver  model.  Equation  2.2, 
the  crossambiguity  function  Is  derived  In  Appendix  A  as 


I  10  1  < 

X(T  ,«  )  -  — 
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slntuYfl2 


(A— 9) 


where  t1  and  »'  are  defined  In  Appendix  A  for  zero  Doppler,  Equation 
A-9  reduces  to  simply  T/Ts,  a  constant. 


The  slgnal-to-nolse  ratio  may  now  be  written  In  terms  of  the 


scattering  and  ambiguity  functions  as. 


(2.21) 


Substituting  Equation  2.16  Into  2.21  yields 
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Equation  2.23  predicts  the  mean  signal  to  mean  noise  ratio  as  a 
function  of  transmit  and  processing  pulse  lengths.  It  may  be 
evaluated  numerically  by  knowing  the  scattering  strengths  pE  and 
E.  for  a  particular  scatterer. 


CHAPTER  3 


SIMULATION  OF  TARGET  ANO  RECEIVER  MODELS 


3.1  Simulation  of  the  Composite  Scattering  Function 


3.1.1  Simulation  of  a  Continuous  Scatterer.  The  model  for  a 

range  spread  scatterer  from  VanTrees  £16]  Is  used  to  simulate  the 

continuous  scatterer.  The  continuous  scatterer  Is  modeled  as  a  line 

of  equl-spaced  point  scatterers,  each  with  a  random  amplitude  drawn 

from  Independent  zero-mean,  complex,  Gaussian  processes  with  variance 
2 

a  .  A  complex  6auss1an  process  Is  defined  as 

?<»»)  -  rR(n)  ♦  j  ?j(n)  (3.1) 


where  rR  and  r^  are  sample  functions  from  real  Gaussian 
variables.  It  Is  assumed  that  ?R  and  r^  are  uncorrelated,  so  that 
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^(njFj 


0,  for  1  *  J  , 


(3.2) 


The  continuous  scatterer  Is  approximated  by  1000  point 
scatterers,  each  of  variance  <*c$,  and  separated  by  l/f$  *  .01 
ms,  where  f  Is  the  digital  sampling  frequency  In  Hertz.  The 
variance  of  each  point  In  the  continuous  scatterer  Is  related  to  the 
continuous  scattering  density,  p£,  by 
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3.1.2  Simulation  of  Point  Scatterers.  The  model  for  the  point 
scattering  component  Is  the  slowly  fluctuating  point  target,  defined 
by  VanTrees  [17].  The  simulation  model  Is  specified  by  fixing  the 
location  of  the  point  scatterers  along  the  length  of  the  continuous 
scatterer.  For  the  purposes  of  this  thesis,  the  point  scatterers  are 
fixed  at  sample  locations  300  and  700,  on  the  1000  sample  continuous 
scatterer.  It  Is  assumed  the  amplitudes  of  the  point  scatterers  are 
sample  functions  of  zero  mean,  uncorrelated  Gaussian  process,  with 
equal  variance,  ojjj,  where  ajjj  >  a ‘  . 

To  calculate  the  energy  reflected  from  the  point  scatterers,  the 
length  of  the  transmit  pulse  Is  restricted  to 


T  <  700-300  -  40  ms  (3.4) 

5  fs 

where  300  and  700  are  the  sample  numbers  corresponding  to  the  location 
of  the  two  point  scatterers.  Equation  3.4  Implies  the  transmit  pulse 
length  Is  short  enough  to  resolve  the  point  reflections,  though  not 
the  scatterers  In  the  continuous  scatterer,  and  the  average  energy  In 
N  point  scatterer  reflections  Is  then  N  times  the  energy  of  a  single 
point  reflection. 


3.2  Generation  of  Simulated  and  Predicted  SNR 

3.2.1  Simulation  of  Reflected  Signal.  To  generate  the  simulated 
signal  reflection  from  the  composite  scatterer,  many  realizations  of 
the  random  scattering  process,  whose  mean  reflection  characteristics 
are  described  by  the  scattering  function,  are  required.  Each  sample 
of  this  random  process  Is  obtained  from  a  Gaussian  random  number 
generator,  described  In  Appendix  B,  assuming  zero  mean,  and  a  standard 
deviation  of  ac$,  or  aN,  depending  on  the  scattering 
component  required.  To  simulate  the  composite  scatterer  In  noise,  the 
outputs  of  three  Independent  random  number  sequences,  one  for  each 
random  component,  are  added  to  form  a  single  realization  of  the  random 
composite  scatterer.  The  envelope  of  the  transmit  signal  Is  the 
convolved  with  this  random  sequence  to  produce  a  representation  of  the 
signal  reflected  from  the  random,  range-spread,  composite  scatterer. 
This  signal  Is  passed  through  the  correlation  receiver,  and  the 
average  of  many  receiver  outputs  compared  to  the  predicted  output. 

The  convolution  Is  accomplished  by  Fourier  transforming  the  transmit 
and  scattering  sequences,  multiplying  them  together,  and  performing  an 
Inverse  Fourier  Transform.  The  size  of  the  FFT  Is  determined  by 


NFFT  -  2  [L0G10(NT_+NTAR)/A10G(2) ]  +  1 


(3.5) 


where  NT^  -  T  «f  an<*  NTAR  “  1000.  The  result  Is  a  sequence 
1099  samples  long,  with  samples  300  to  400  and-700  to  800  representing 
the  point  reflector  returns.  Table  1  lists  the  simulation  parameters 
that  are  constant  for  all  runs.  Including  the  durations  of  the 
transmit  signal  and  continuous  scatterer. 

Noise  Is  assumed  to  be  zero  mean,  white  and  Gaussian.  It  Is 
simulated  Identically  to  the  continuous  scattering  sequence  discussed 
In  Section  3.1.1,  with  the  noise  spectral  density  equal  to  the 
variance  of  each  point  In  the  random  noise  sequence. 

The  simulation  of  sequences  of  random  numbers,  both  for  the 
scattering  process  and  noise  1$  discussed  In  Appeddlx  B. 


3.2.2  Generation  of  Simulated  and  Predicted  SNR.  To  calculate 
the  simulated  SNR,  each  realization  of  the  reflected  signal  Is  passed 
through  the  digital  equivalent  of  the  receiver  shown  In  Figure  1,  for 
a  given  number  of  samples  NT,  where  NT  *  T*f$.  J  realizations  of 
the  signal,  and  J  realizations  of  noise  are  passed  through  the 
receiver  separately,  with  the  mean  slgnal-to-nolse  ratio  computed  by 


Table  1  -  Constant  Simulation  Scattering  and  Receiver  Model  Parameters. 
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-  Continuous  scatterer  length  *  100  ms 

-  Transmit  pulse  length  -  10  ms 

-  Sampling  rate  -  10000  Hz 

-  Number  of  highlight  scatterers  *  2 

-  Location  of  highlights  *  30  and  70  ms 

-  Minimum  Receiver  Integration  Time  -  10  ms 


b 

V 


/  4 


where  S 


Is  the  K—  receiver  output  for  signal  only  and 


/th 


N^  ‘  Is  the  K—  output  for  noise  only. 


The  predicted  SNR  Is  evaluated  by  relating  the  terms  In  the 
SNR  prediction  equation.  Equation  2.23,  to  the  variance  of  each  of  the 
scattering  and  noise  variances,  or 


SNR 


v[ 


2.2—  s  N  2 
*CS  fs  Ts  +  T  Nhi  °HI 


] 


(3.7) 


where  N  Is  the  number  of  samples  of  the  point  scatterer  returned  In 


the  Interval  [0,T]  and  0  <  N/NTS  <  2.  Is  the  variance  of  the 


zero  mean,  Gaussian,  noise  process. 


As  noted  earlier,  the  range  of  T/Ts  values  discussed  In  the 
thesis  are  between  1  and  11.  For  simplicity.  It  is  assumed  that  any 
bulk  delay  In  the  reflected  signal  Is  known,  so  that  the  receiver 
begins  processing  the  reflected  signal  at  the  onset  of  the  continuous 
scatterer  return,  T/T$  -  0.  The  first  comparison  between  simulation 
and  Equation  3.7  Is  made  at  T/Ts  -  1,  after  an  entire  transmit  pulse 
duration  has  been  processed  through  the  correlation  receiver.  This 
method  of  fixing  the  processing  start  time  is  followed  throughout  the 


thesis,  whether  the  continuous  scatterer  return  Is  present  or  not.  In 
particular,  the  examples  shown  In  Section  4.3,  point  scatterers  in 


noise,  begin  at  T/Ts  «  4,  one  transmit  pulse  length  after  the  first 
point  reflector  return  at  T/Ts  -  3,  with  the  receiver  output  from 
T/Ts  *  1  to  T/Ts  «  2.5  representing  noise  only. 

In  addition,  the  energy  in  the  signal  reflected  from  the 
continuous  scatterer  Is  not  constant  with  Increasing  Integration  time, 
as  Implied  by  In  Equation  3.3,  and  the  relationship  of  the  continuous 
scattering  variance  to  the  continuous  scattering  density  must  be 
modified  to  account  for  the  time  spread  of  the  reflected  signal.  The 
received  signal  Is  the  convolution  of  the  transmit  signal,  say  N 
points  long,  and  the  target,  H  points  long,  where  N  >  N.  The  resulting 
sequence  Is  M+N-l  long.  In  Figure  2,  note  that  the  received  signal 
duration  is  roughly  0.11  sec,  for  a  scattering  function  0.1  sec  long 
and  a  0.01  sec  transmit  pulse.  The  variance  of  the  signal,  assuming 
uncorrelated  point  scatterers  Is 


The  variance  of  the  first  N  samples  of  the  reflected  signal,  for 

3  2 

a  transmit  pulse  length  of  N  samples.  Is  not  N  •  ac$  but, 

<’E06E-("3-2("t,»  °cs2  (3'10 

Equations  3.9  and  3.10  are  then  used  as  corrections  for  calculations 
of  the  continuous  scattering  energy  In  Equation  3.7. 


CHAPTER  4 


COMPARISON  OF  PREDICTED  AND  SIMULATED  SNR 


Introduction 


4.1,1  Introduction.  In  this  chapter,  the  comparison  of  the 
predicted  and  simulated  SNR  Is  presented.  Fourteen  combinations  of 
scattering  strength  parameters  are  presented.  These  cases  cover  small 
and  large  point  scatterer  reflectivities,  and  both  high  and  low 
slgnal-to-nolse.  The  standard  deviation  of  the  scattering  strength 
parameters  are  listed  In  Table  2,  along  with  some  statistics  of  the 
differences  between  predicted  and  simulated  SNR.  The  constants  for 
all  the  simulation  trials  are  transmit  pulse  length,  T  ,  continuous 
scatterer  duration,  point  scatterer  position  along  the  continuous 
scatterer  and  the  number  of  Independent  realizations  of  the  model 
averaged  to  compare  with  the  SNR  predicted  from  Equation  3.8. 


Table  2  lists  the  scattering  strength  parameters  and  some 
statistics  of  the  differences  between  predicted  and  average  simulation 
SNR,  for  the  fourteen  scattering  strength  parameters  examined. 

Columns  2-4  list  the  standard  deviation  of  the  scattering  and  noise 
components  of  Equation  3.7.  Columns  5-7  list  the  maximum  difference 
between  simulation  and  prediction,  the  geometric  mean  difference 
between  prediction  and  simulation,  and  the  standard  deviation  of  the 
difference,  respectively,  for  each  set  of  scattering  strength 
parameters. 


TABLE  2  -  Scattering  Strength  Parameters  and  Statistics  of  Differences 
between  Simulated  and  Predicted  and  Simulated  Signal -to-Noise  Ratios 
for  14  Combinations  of  Scattering  Strength  Parameters. 


STANDARO  DEVIATION  OF 

SCATTERING  DIFFERENCES  BETWEEN 

PARAMETERS  MODEL  AND  THEORY 

(Average  of  all  values 
of  T/Ts.) 


CASE  NO. 

°CS 

°HI 

°N 

MAX 

MEAN 

ST.  DEV. 

DB 

DB 

DB 

1 

1. 

0. 

1. 

1.1 

0.4 

0.3 

2 

0. 

1. 

1. 

1.0 

0.4 

0.3 

3 

100.. 

0. 

1. 

1.1 

0.5 

0.3 

4 

0. 

100. 

1. 

1.2 

0.5 

0.3 

5 

1. 

0. 

100. 

1.6 

0.6 

0.4 

6 

0. 

1. 

100. 

1.7 

0.4 

0.4 

7 

1. 

10. 

1. 

1.2 

0.5 

0.3 

8 

1. 

100. 

1. 

1.3 

0.5 

0.3 

9 

1. 

100. 

100. 

1.4 

0.6 

0.3 

10 

1. 

1. 

1. 

1.5 

0.5 

0.4 

11 

1. 

0. 

1000. 

0.9 

0.3 

0.3 

12 

1. 

100. 

1000. 

2.2 

0.4 

0.5 

13 

1. 

0. 

9500. 

1.4 

0.6 

0.4 

14 

1 . 

100. 

9500. 

1.1 

0.4 

0.4 

AVERAGE 

1.3 

0.5 

0.3 

The  mean  difference  between  theory  and  Equation  3.7  Is  defined  as 
the  geometric  average  of  differences  at  each  value  of  T/T$  examined. 

In  decibels,  for  a  set  of  scattering  strength  parameters,  or, 

i  N 

f(d8)  ”  n  L  E1  (dB)  (4J) 

1-1 

where  Ej  Is  the  absolute  value  of  the  difference  of  the  predicted 
SNR  and  the  simulated  SNR,  In  dB,  for  a  particular  value  of  T/T$. 
Similarly,  the  standard  deviation  of  differences  Is  computed  by 


scattering  strength  parameters  a  ,  oUT,  ou.  Equation  3.7  Is 

CS  HI  N 

used  to  generate  the  theoretical  SNR  prediction  for  receiver 
Integration  times  from  equal  to  the  transmit  pulse  length,  T/T$  - 
1.0,  to  T/T$  -  11,  which  encompasses  the  entire  scattered  signal. 

The  prediction  equation  Is  evaluated  at  Increments  of  one-half  the 
transmit  pulse  length  between  these  two  values. 

The  simulation  SNR  Is  also  generated  for  each  set  of  scattering 
strength  parameters,  and  the  21  values  of  T/T$  between  1.0  and  11. 
100  realizations  of  the  random  model,  which  are  assumed  to  be 
Independent  since  each  starts  with  a  unique  seed  value,  are  linearly 


averaged.  Ten  times  the  logarithm  of  the  average  value  Is  then 
compared  with  the  predicted  value. 


It  Is  Important  to  note  that  the  statistics  of  the  differences 
between  theory  and  model  listed  In  Table  2  are  calculated  using  only 
19  of  the  receiver  Integration  times  between  T/T$  ■  1.0  and  11.  Due 
to  the  point  scatterer  return  at  T/T$  ■  3.0,  the  differences  between 
model  and  prediction  equation  at  T/T$  «  3.0  and  3.5  are  not 
Included.  At  these  Integration  times,  the  received  signal  Includes  a 
statistically  small  number  of  samples  from  the  large  variance, 
highlight  return,  and  therefore  displays  large  fluctuations,  even 
after  averaging  100  simulation  realizations.  Not  until  a  complete 
transmit  pulse  Is  processed  by  the  receiver.  Is  a  comparison  between 
the  predicted  result,  and  the  scattering  model  output,  considered 
valid. 

4.1.3  Presentation  of  Comparison  Data.  In  addition  to  the 
results  listed  In  Table  2,  two  types  of  graphical  comparisons  are 
used.  An  example  of  the  first  plot  type  Is  shown  In  Figure  3,  with 
SNR  on  the  ordinate,  and  T/Ts  on  the  absclcca.  The  theoretical  SNR 
prediction  Is  shown  as  circles,  and  the  model  average  SNR  depicted  by 
squares.  Note  the  data  points  at  T/Ts  -  3.0  and  3.5  are  shown  In 
this  figure,  though  they  are  not  Included  In  the  statistical 
comparisons  listed  In  Table  2. 


The  second  type  of  graph  Is  shown  In  Figure  4.  As  In  Figure  3, 
the  predicted  values  are  shown  as  circles,  and  the  simulation  data  as 
squares.  The  absclcca  axis  Is  the  ratio  of  transmit  pulse  length  to 


N 


receiver  Integration  time,  T/T$.  The  ordinate  Is  SNR,  normalized  by 
the  predicted  SNR.  Thus  the  predicted  values  appear  to  lie  on  a 
horizontal  line  at  0  dB*  The  vertical  bars  around  the  normalized 
simulation  data  depict  the  90  percent  confidence  Interval  to  be 
discussed  In  Section  4.2. 
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4.2  The  Continuous  Scatterer  In  Noise 

4.2.1  Continuous  Scatterer  In  Noise.  The  simplest  realization  of 
the  model  Is  the  continuous  scatterer  only  In  the  presence  of  low 
background  noise.  Figure  3  Is  a  SNR  versus  T/Ts  plot,  for  Case  1,  a 
high  SNR  combination  of  scattering  strength  parameters  (see  Table  2). 
The  maximum  predicted  SNR  occurs  at  a  value  of  T/T$  -  10.5.  The 
simulation  SNR  values  generally  follow  the  prediction  curve,  with  a 
maximum  difference  between  the  two  curves  of  1.1  dB,  and  a  mean 
difference  of  0.4  dB. 


E 


k 


The  prediction  curve  has  maximum  slope  for  T/T$  values  less 


than  4.0.  The  slope  of  the  SNR  versus  T/Ts  curve  Is  Important, 


since  It  Indicates  the  region  of  receiver  Integration  times  that  yield 
the  greatest  Increase  In  SNR,  for  the  least  Increase  In  Integration 
time.  This  behavior  Is  Intuitively  correct,  since  Increasing  the 


Integration  time  from  T/T  -  1.0  to  2.0,  doubles  the  amount  of 
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signal  energy  processed  by  the  receiver.  Conversely,  Increasing  the 


listen  "window"  from  T/Ts  ■  10.0  to  11.0  adds  only  ten  percent  more 
signal  energy. 

4.2.2  Distribution  of  Simulated  SNR.  The  ambiguity  -  scattering 
function  analysis  presented  In  Chapter  2  allows  calculation  of  the 
mean  SNR  for  a  given  set  of  scattering  parameters  and  pulse  lengths, 
but  contains  no  Information  concerning  the  distribution  of  SNR  about 
the  mean.  Using  VanTrees*  models  for  noise  and  the  continuous 
scatterer.  It  Is  possible  to  derive  the  distribution  of  simulated  SNR. 

The  essence  of  the  derivation  Is  that  both  received  signal  and 
noise  are  modeled  as  Independent,  zero  mean,  6auss1an  processes.  The 
quantity  of  Interest,  SNR,  Is  the  ratio  of  received  signal  squared  and 
noise  squared.  The  SNR  Is  then  the  ratio  of  two  chi-squared 
variables,  and  the  distribution  of  this  ratio  Is  the  F  distribution. 
The  Interested  reader  Is  referred  to  Appendix  C  for  the  complete 
derivation. 

Knowing  the  distribution  of  SNR,  the  confidence  limits  of  the 
simulated  data  may  readily  be  computed.  Figure  4  Is  a  plot  of 
normalized  SNR  versus  T/Ts  for  the  case  shown  In  Figure  3.  As  noted 
previously,  both  the  simulated  and  predicted  SNR  values  are  normalized 
by  the  predicted  SNR,  which  transform  the  prediction  curve  In  Figure  3 
to  the  horizontal  line  In  Figure  4.  For  90  percent  confidence,  the 
simulated  values  must  lie  within  +  0.75  dB  of  the  predicted  value. 


This  range  Is  depicted  In  Figure  4  by  the  vertical  bars  through  each 
of  the  normalized  simulation  data.  Note  lit  Figure  4  that  three  of  the 
twenty-two  normalized  simulation  values  are  not  within  the  confidence 
bounds,  roughly  the  number  expected  for  90  percent  confidence. 

Also  It  Is  worth  noting  that  the  confidence  Intervals  for  the 
matched  filter  condition  are  not  strictly  correct.  For  T/T$  ■  1,  the 
received  signal  builds  from  one  sample  at  t  *  0  to  the  sum  of  N 
scatterers  at  t  •  T  .  The  sequence  Is  triangular,  with  each  sample 
Increasing  In  degrees  of  freedom  with  time.  The  total  number  of 
degrees  of  freedom  In  this  sequence  are  not  the  same  as  the  rest  of 
the  received  signal,  which  results  In  a  slightly  wider  confidence 
Interval  than  the  ±  0.75  dB  for  values  of  T/T$  w  1. 

4.2.3  Comparison  as  a  Function  of  SNR.  The  comparison  between 
simulation  and  prediction  appears  to  be  unaffected  by  the  absolute 
SNR.  In  addition  to  the  case  discussed  above,  three  other 
combinations  of  continuous  scattering  strength  and  noise  were 
examined,  with  the  predicted  SNR  between  -20  to  +95  dB.  Representative 
examples  are  shown  In  Figures  5  and  6. 


Note  Figures  5  and  6  Indicate  only  small  differences  between 
model  and  theory  for  the  case  of  equal  continuous  scattering  strength 
and  noise  spectral  density,  l.e.,  0  dB.  The  largest  difference 
between  prediction  and  simulation  Is  1.3  dB,  and  as  can  be  seen  In 


I 


I  FIGURE  5.  Comparison  of  predicted  and  simulated  SNR  for 

Case  11 ;  ocs  -  1 .  oHI  ■  °-  °N  *  100°-  Predicted 
►  values  are  shown  as  circles,  simulated  data  shown 
*  as  squares. 


NORMALIZED  SNR  (dB) 


Figure  6  only  two  of  the  data  points,  at  T/Ts  *  1.5  and  3.0,  fall 
outside  the  90  percent  confidence  Interval.  Similar  results  were 
observed  for  for  both  very  low  SNR,  -20  dB,  and  very  high  SNR,  +95  dB. 

4.3  Point  Reflectors  In  Noise 

The  other  components  of  the  composite  scatterer  are  the  two  point 
reflectors.  Figure  7  Is  a  plot  of  SNR  versus  T/Ts,  for  two 
moderately  large  point  scatterer  returns  In  noise.  Since  the  point 
scatterer  return  does  not  begin  until  T/Ts  -  3.0,  comparisons 
between  the  model  and  theory  are  not  considered  valid  until 
T/T$  -  4.0.  Note  that  absclcca  values  In  Figures  7-10  begin  at  4.0 
for  this  reason. 

Figure  7  Is  a  SNR  versus  T/Ts  plot  for  the  two  highlight 
reflectors  In  noise  with  a  mean  SNR  of  roughly  35  dB.  As  before,  the 
predicted  values  are  shown  as  circles,  and  the  simulation  values  as 
squares  or  boxes.  The  largest  difference  between  the  two  curves  Is 
1.0  dB,  with  the  average  difference  0.4  dB,  Indicating  that  the  model 
accurately  follows  the  predicted  values.  Note  that  both  curves  are  a 
maximum  at  T/Ts  *  4.0  and  8.0,  corresponding  to  receiver  integration 
times  that  receive  the  entire  point  scatterer  reflection  and  minimum 
noise  energy.  Between  the  first  and  second  point  reflections,  and 
after  the  second  point  reflection,  the  SNR  drops  off  due  to  the 
addition  of  noise,  and  the  lack  of  any  additional  signal  energy. 


FIGURE  7.  Comparison  of  predicted  and  simulated  SNR  for 
Case  2;  ocs  *  0,  *  1 ,  oN  -  1 .  Predicted 

values  are  shown  as  circles,  simulated  data  shown 
as  squares. 


The  good  agreement  between  theory  and  model  is  also  observed  in 


the  normalized  SNR  versus  T/T$  plot.  Figure  8.  Three  of  16 
simulation  values  fall  outside  the  confidence  Interval.  No  systematic 
or  periodic  differences  between  the  two  curves  is  apparent. 

A  low  SNR  case  is  shown  In  Figures  9  and  10.  Note  the  rapid 
decrease  In  SNR  between  the  point  reflector  returns  when  no  continuous 
scattering  component  Is  present.  Again,  the  simulated  values  of  SNR 
closely  follow  the  predicted  values,  and  the  normalized  data.  Figure 
10,  shows  only  a  single  point  out  of  16  outside  the  90  percent 
confidence  Interval. 

4.4  The  Composite  Scatterer  In  Noise 

4,4.1  Composite  Scatter  In  Noise.  A  complete  realization  of  the 

model  includes  reflections  from  two  point  reflectors,  of  variance 
2 

oHI,to  the  continuous  scatterer  return.  The  locations  of  the 
point  scatterers  are  fixed  to  T/T$  values  of  3.0  and  7.0.  Two 
examples  of  the  envelope  of  the  signal  reflected  from  the  composite 
scatterer  are  shown  In  Figures  2  and  11.  The  point  scattering 
strength  In  Figure  2  Is  10  times  the  point  scatterer  strength  in 
Figure  11,  but  In  both  cases,  the  variance  of  the  point  scatterers  Is 
large  compared  to  the  continuous  scatterer  and  noise. 


FIGURE  9.  Comparison  of  predicted  and  simulated  SNR  for 
Case  6;  ocs  «  0,  oHI  «  1,  -  TOO.  Predicted 

values  are  shown  as  circles,  simulated  data  shown 
as  squares. 


FIGURE  10.  Comparison  of  normalized  SNR  for  Case  6;  ocs  *  0, 
°HI  s  1 •  ou  *  100.  Predicted  values  are  shown  as 
circles,  simulated  data  shown  as  squares. 


TIME  (sec) 


FIGURE  11.  One  realization  of  the  signal  reflected  from  the 
composite  scatterer  for  large  point  scatterer 
reflectivity. 


The  SNR  versus  T/T$  plot  for  a  high  SNR,  large  point  scatterer 
reflectivity  case  is  shown  In  Figure  12.  For  values  of  T/T$  less 
than  3.0,  the  signal  energy  Is  due  soley  to  the  continuous  scatterer. 
With  the  onset  of  the  first  point  scatterer  reflection,  the  SNR 
rapidly  Increases,  with  the  prediction  curve  reaching  a  local  maximum 
at  T/T$  «  4.0.  Between  the  point  scatterer  returns,  the  SNR  falls 
off  3  dB,  and  reaches  a  second  maximum  with  the  arrival  of  the  second 
point  reflection.  Differences  between  simulation  and  prediction  for 
this  example  are  roughly  equal  to  the  differences  seen  In  the  previous 
sections.  Again  neglecting  data  at  T/T^  ■  3.0  and  3.5,  the  average 
difference  between  the  two  curves  Is  0.5  dB,  with  a  maximum  difference 
of  1.3  dB,  at  a  value  of  T/Ts  *  10.0. 

The  normalized  SNR  plot  for  this  example  Is  shown  In  Figure  13. 
The  majority  of  the  normalized  simulation  values  fall  within  the  90 
percent  confidence  Interval.  Only  the  data  at  T/T$  -  6.0  and  10.0 
are  outsides  these  limits.  It  Is  Important  to  note  that  the  composite 
scatterer  data,  as  well  as  the  component  scatterers  Individually,  does 
not  display  any  apparent  pattern  In  the  differences  between  simulation 
and  prediction. 

Figure  14  Is  a  SNR  versus  T/Ts  plot  of  composite  scatterer  data 
for  point  scatterer  reflectivities  1/10  the  strength  of  the  example 
shown  In  Figure  12.  The  predicted  SNR  values  for  T/T$  less  than  3.0 
are  the  same  In  both  Figures  12  and  14,  though  plotted  on  different 
scales.  The  onset  of  the  first  point  scatterer  reflection  In  Figure 


FIGURE  12.  Comparison  of  predicted  and  simulated  SNR  for 
Case  8;  ocs  *  1,  *  100,  oN  *  1 .  Predicted 

values  are  shown  as  circles,  simulated  data  shown 
as  squares. 
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FIGURE  14.  Comparison  of  predicted  and  simulated  SNR  for 
Case  7;  ocs  ■  1,  -10,  *  1 .  Predicted 

values  are  shown  as  circles,  simulated  data  shown 
as  squares. 


14,  smaller  reflectivity.  Is  less  dramatic  than  In  Figure  12,  larger 
reflectivity,  with  the  predicted  SNR  Increasing  only  2  dB,  over  the 
continuous  scatterer  level.  The  prediction  curve  In  Figure  14  does 
decrease  slightly  between  the  point  scatterer  returns,  as  expected, 
and  jumps  up  at  the  start  of  the  second  point  scatterer  return.  Note 
that  the  different  scales  between  Figures  12  and  14  give  the  - 
appearance  that  the  simulation  values  are  In  Figure  14  do  not  follow 
the  prediction  curve  as  closely  as  the  example  shown  In  Figure  12. 
However,  the  differences  between  theory  and  model  are  roughly  equal  to 
the  previous  example,  as  shown  In  Table  2. 

The  normalized  SNR  versus  T/T$  plot  for  this  set  of  scattering 
parameters  Is  shown  In  Figure  15.  Four  of  the  19  valid  simulation 
normalized  SNR  values  are  outside  the  90  percent  confidence  limits, 
but  the  simulation  values  are  more  closely  aligned  with  the  normalized 
prediction  values  than  the  data  shown  In  Figure  13. 

4.4.2  Effects  of  Simulation  Sample  Size.  All  of  the  simulation 
data  presented  up  to  this  point  was  obtained  by  averaging  100 
simulation  realizations.  The  differences  between  model  and  theory  may 
be  reduced  by  Increasing  the  number  of  averages  used  In  the 
comparisons.  Figure  16  Is  a  SNR  versus  T/Ts  plot  for  the  set  of 
scattering  parameters  shown  In  Figure  12,  but  obtained  using  the 


FIGURE  16.  Comparison  of  predicted  and  simulated  SNR  for 
Case  8;  ocs  -1,  ouj  *  100,  om  -  1 .  Simulated 
SNR  values  obtained  by  averaging  1000  realizations 
of  the  scattering  function.  Predicted  values  are 
shown  as  circles,  simulated  data  shown  as  squares. 


average  of  1000  simulation  trials  Instead  of  100.  The  prediction 
curve  values  are  Identical  between  this  figure  and  Figure  12,  as 
expected.  However,  the  simulation  values  lie  much  more  closely  to  the 
prediction  curve  for  the  1000  trial  average  than  for  only  100  trials. 
The  mean  difference  between  the  theory  and  model  values  In  Figure  16, 
Is  only  0.2  d8,  compared  to  0.5  dB,  for  Figure  12. 

The  excellent  agreement  between  the  prediction  equation  and  model 
Is  also  observed  In  the  normalized  SNR  versus  T/Ts  plot.  Figure  17. 
Only  the  data  point  at  T/T$  -  3.5  lies  outside  the  90  percent 


confidence  limits. 
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CHAPTER  5 


RECEIVER  PERFORMANCE  FOR  CORRELATED  CONTINUOUS  SCATTER 


5.1  Modifications  to  VanTrees1  Model 

5.1 .1  Introduction.  The  model  of  the  range  spread  scatterer, 
formulated  by  VanTrees,  was  simulated  as  a  set  of  equl-spaced  point 
reflections.  The  reflectivity  of  these  point  scatterers  Is  modeled  as 
sample  functions  from  uncorrelated  Gaussian  processes.  In  this 
chapter  the  model  Is  modified  to  allow  varying  degrees  of 
Interdependence,  or  correlation,  between  the  Individual  point 
reflections  that  comprise  the  continuous  scattering  return. 

This  modification  violates  the  conditions  under  which  the 
scattering  function  Is  defined  In  Chapter  2.  Indeed,  for  correlated 
scattering,  the  scattering  function  must  be  replaced  with  the  more 
general  time-frequency  correlation  function  [18],  By  comparing 
simulated  SNR,  produced  using  the  correlated  continuous  scatterer,  and 
the  predicted  SNR,  which  assumes  uncorrelated  scattering,  the  effects 
of  correlated  scattering  on  the  validity  of  the  prediction  may  be 
examined. 


5.1.2  Simulation  of  the  Correlated  Continuous  Scatterer.  The 


continuous  scattering  components  may  then  be  achieved  by  low  or  band 
pass  filtering  the  continuous  scattering  sequence.  The  degree  of 
correlation  will  be  Inversely  proportional  to  the  filter  bandwidth. 

Eight  filter  bandwldths  are  used  to  generate  varying  degrees  of 
correlation.  In  normalized  frequency,  transmit  frequency  divided  by 
the  sampling  frequency,  the  normalized  filter  bandwldths  range  from 
0.49,  almost  uncorrelated,  to  0.01,  almost  perfectly  correlated.  The 
actual  shape  of  the  filter  response  Is  relatively  unimportant,  as  Is 
the  difference  between  low  and  band  pass.  The  actual  filters  used  are 
based  on  a  low  pass,  cosine  tapered,  frequency  sampling  filter  [20]. 
Figures  18  and  19  are  examples  of  the  filter  frequency  response  for 
normalized  cutoff  frequencies  of  0.49  and  0.01. 

To  produce  the  correlated  scattering  sequence,  the  uncorrelated 
sequence  Is  convolved  with  the  filter  Impulse  response.  The  filter 
output  Is  then  normalized  by  the  ratio  of  the  Input  sequence  energy  to 
the  energy  In  the  filtered  sequence,  producing  a  final,  correlated 
scatterer,  with  the  same  total  energy  as  the  original  uncorrelated 
sequence.  This  normalization  Insures  that  changes  In  the  receiver 
output  SNR  are  due  solely  to  the  correlation  among  the  continuous 
scattering  components. 

Figures  20-24  are  plots  of  an  unfiltered  sequence  and  the  output 
of  four  filters.  In  decreasing  bandwidth.  Figure  19,  generated  using 
a  filter  bandwidth  of  0.49  is  Indistinguishable  from  the  unfiltered 
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FIGURE  18.  Frequency  response  of  filter  used  to  generate  the 

correlated  scattering  sequence.  The  filter  bandwidth 
is  0.49. 
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FIGURE  23.  Correlated  scattering  sequence  for 
of  0.1.  Sequence  energy  *  65.95  dB 


FIGURE  24.  Correlated  scattering  sequence  for  a  filter  bandwidth 
of  0.01.  Sequence  energy  *  65.95  dB. 


data.  As  the  bandwidth  shrinks,  the  sequence  appears  more  sinusoidal 
than  random,  as  expected.  Note  the  energy  of  each  sequence  In 
Figures  20-24  Is  Identical,  65-95  dB.  To  generate  the  correlated 
composite  scatterer,  the  two  point  reflectors  are  then  added  to  the 
filtered  continuous  scattering  sequence  at  samples  300  and  700,  as 
discussed  previously. 


5.2  Correlated  Continuous  Scatterer  Results 

5.2.1  SNR  Comparison  versus  Bandwidth.  Table  3  lists  the 
comparison  statistics  for  the  correlated  scatterer  simulation  and 
predicted  uncorrelated  SNR.  The  related  figures  are  25-29.  As  with 
the  uncorrelated  continuous  scatterer  data,  comparisons  between  the 
simulated  and  predicted  SNR  values,  do  not  Include  data  at  T/T$  = 
3.0  and  3.5. 

In  Figure  25,  the  filter  bandwidth  Is  large,  and  the  comparison 
between  simulation  and  prediction  Is  quite  good.  Indeed,  the 
statistical  comparison  for  this  data  set  Is  as  good  as  any 
uncorrelated  data  comparison.  The  quality  of  the  comparison  Is  also 
seen  In  the  normalized  comparison.  Figure  26.  Note  the  majority  of 
the  prediction  values  fall  within  the  90  percent  confidence  Interval; 
only  the  values  at  T/T$  *  5.5  and  7.0  are  outside  this  range.  The 
average  difference  between  simulation  and  prediction  Is  1.1  dB, 
comparable  to  all  the  uncorrelated  scattering  data. 


TABLE  3  -  Scattering  Strength  Parameters  and  Statistics  of  Differences 
between  Simulated  and  Predicted  and  Simulated  Slgnal-to-Nolse  Ratios 
Assuming  Correlated  Continuous  Scattering. 


STANDARD  DEVIATION  OF 

SCATTERING  MODEL  DIFFERENCES  BETWEEN 

PARAMETERS  MODEL  AND  THEORY 

(Average  of  all  values 
of  T/Ts.) 


°N 


CASE  NO.  BW 


«es 


®HT 


MAX  MEAN  ST.  DEV. 


FIGURE  25.  Comparison  of  uncorrelated  continuous  scattering  SNR 
prediction,  and  correlated  scattering  simulation  for 
Case  15;  oCe  ■  1,  or  *  1 ,  filter  bandwidth 

0.49.  Predicted  values  are  shown  as  circles, 
simulated  data  shown  as  squares. 
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FIGURE  26.  Comparison  of  normalized  uncorrelated  continuous 

scattering  SNR  prediction,  and  correlated  scattering 
simulation  for  Case  15;  ocs  «  1,  *  0,  oN  *  1 , 

filter  bandwidth  =  0.49.  Predicted  values  are  shown  as 
circles,  simulated  data  shown  as  squares. 


FIGURE  28.  Comparison  of  uncorrelated  continuous  scattering  SNR 
prediction,  and  correlated  scattering  simulation  for 
Case  17;  oCc  -  1,  oHt  *  0,  oN  *=  1,  filter  bandwidth 
*  0.1.  Predicted  values  are  shown  as  circles, 
simulated  data  shown  as  squares. 


FIGURE  29.  Comparison  of  uncorrelated  continuous  scattering  SNR 
prediction,  and  correlated  scattering  simulation  for 
Case  18;  ocs  =  1,  <jhi  *  0,  or  =  1 ,  filter  bandwidth 
«  0.01.  Predicted  values  are  shown  as  circles, 
simulated  data  shown  as  squares. 
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5.2.2  Corrections  for  Correlated  Scattering.  The  SNR  comparisons 
In  Figures  25-29  suggest  a  correction  factor  might  be  applied  to  the 
uncorrelated  prediction  to  account  for  the  presence  of  correlated 
scattering.  Two  different  types  of  curve  fits  are  employed.  The 
first  Is  suggested  by  the  simulation  data,  and  assumes  the  correction 
equation  is  of  the  form, 

SNR  Gain  (dB)  =  ea(0,5~f)  (5.1) 

where  SNR  Gain  Is  the  mean  difference  between  the  predicted  SNR  values 
and  correlated  simulation  SNR;  f  Is  the  normalized  filter  bandwidth, 
and  it  Is  assumed  the  SNR  Gain  at  0.5  Is  zero.  Several  values  of  the 
coefficient  a  were  plotted  to  obtain  a  empirical  best  fit.  Figure  30 
is  a  plot  of  the  SNR  Gain  versus  filter  bandwidth  for  the  eight  cases 
examined.  Also  plotted  are  three  realizations  of  Equation  5.1,  with 
values  of  a  *  5.6,  5.8  and  5.9.  All  three  curves  provide  a  reasonable 
fit  for  frequencies  greater  than  0.25.  At  frequencies  less  than  0.25, 
all  the  the  fitted  curves  predict  too  large  a  SNR  Gain,  by  1-2  dB. 

For  very  small  filter  bandwldths,  the  gain  indicated  by  the  data 
appears  to  increase  much  faster  than  any  of  the  fits  provided  by 
Equation  5.1. 

A  second  fit  to  the  data  Is  a  simple  polynomial  fit.  For  only  8 
data  points,  the  degree  of  the  fitted  polynomial  is  restricted  to 
between  2  and  4.  Table  4  lists  the  coefficients  computed  by  the  IMSL 
routines  RLFOTH  and  RLDOPH  for  the  three  orders  examined  [21].  A  plot 
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FIGURE  30.  Exponential  fit  to  SNR  gain  versus  filter  bandwidth 
for  a  «  5.6,  5.8  and  5.9.  SNR  gain  obtained  from 
simulation  Is  shown  as  the  solid  line. 
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TABLE  4  -  Coefficients  of  Polynomials  Fitted  to  Continuous  Correlated 
Scatterer  SNR  Sain  versus  Filter  Bandwidth  Data  (See  Figure  33). 


POLYNOMIAL 

DEGREE 

AO 

A1 

A2 

2 

15.7 

-80.9 

105.1 

3 

17.2 

-132.8 

401.4 

4 

18.3 

-194.5 

1025.5 

-411.3 

-2434.8 
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of  the  data  with  the  three  fitted  polynomials  Is  shown  In  Figure  31. 
Clearly,  the  accuracy  of  the  fit  Increases  with  Increasing  polynomial 
degree,  though  each  Increase  In  order  adds  an  additional  Inflection 
point  not  Indicated  by  the  data. 


Figure  32  plots  the  best  mean  square  fit  for  both  models  and  the 
SNR  Gain  data.  The  quality  of  the  polynomial  fit  Is  clearly  better 
than  the  exponential  model,  particularly  In  the  frequency  range  from 
0.2  to  0.05  Hz.  The  fourth  order  polynomial  equation  appears  to  be 
the  most  accurate  fit  to  the  correlated  scattering  gain,  and  might  be 
used  with  the  uncorrelated  scattering  prediction  to  estimate  the 
correlation  of  an  actual  scattering  process. 


5.3  Composite  Correlated  Scatterer  SNR  Comparison 


Four  cases  examine  the  correlated  continuous  scatterer  with  two 
large  amplitude,  point  reflectors.  Table  3  lists  the  simulation  and 
uncorrelated  scattering  prediction  comparison  statistics,  for  constant 
scattering  strength  parameters  and  Increasing  degrees  of  continuous 
scatterer  correlation.  The  related  figures  are  33-38. 


In  Figure  33,  for  a  filter  bandwidth  of  0.49,  the  comparison 
between  simulated  and  predicted  slgnal-to-nolse  ratios  Is  quite  good. 
As  with  the  uncorrelated  scatterer  data,  the  greatest  difference 
between  simulation  and  prediction  occurs  at  T/T$  =  3.5,  an 
Indication  that  the  SNR  Is  dominated  by  the  point  scatterers  returns 


NORMALIZED  BANDWIDTH 


FIGURE  31.  Polynomial  fit  to  SNR  gain  versus  filter  bandwidth 

for  polynomial  orders,  N,  =  2,  3  and  4.  True  SNR  gai 
shown  as  solid  line. 
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FIGURE  32.  Best  fits  for  polynomial  and  exponential  models  to  SNR 
gain  versus  filter  bandwidth.  True  SNR  gain  shown  as 
solid  line. 


FIGURE  33.  Comparison  of  uncorrelated  continuous  scattering  SNR 
prediction,  and  correlated  scattering  simulation  for 
Case  19;  ocs  ■  1,  -  100,  -  1 ,  filter 

bandwidth  *  0.49.  Predicted  values  are  shown  as 
circles,  simulated  data  shown  as  squares. 
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FIGURE  34.  Comparison  of  normalized  uncorrelated  continuous 

scattering  SNR  prediction,  and  correlated  scattering 
simulation  for  Case  19;  ocs  -  1,  o^I  *  °N  *  1 • 
filter  bandwidth  *  0.49.  Predicted  values  are  shown 
as  circles,  simulated  data  shown  as  squares. 


FIGURE  35.  Comparison  of  uncorrelated  continuous  scattering  SNR 
prediction,  and  correlated  scattering  simulation  for 
Case  20;  “1#  «hi  *  100,  ok  *  1,  filter 

bandwidth  *  0.25.  Predicted  values  are  shown  as 
circles,  simulated  data  shown  as  squares. 
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FIGURE  38.  Comparison  of  uncorrelated  continuous  scattering  SNR 
prediction,  and  correlated  scattering  simulation  for 
Case  22;  *  1*  cm  ■  100,  <?#  *  1 ,  filter 

bandwidth  *  0.01.  Predicted  values  are  shown  as 
circles,  simulated  data  shown  as  squares. 


after  T/Ts  =  3.0.  For  this  case,  any  offset  Introduced  by  the 
correlated  continuous  scatterer  Is  smaller  than  the  variation  of  the 
uncorrelated  scatter.  The  normalized  SNR  comparison.  Figure  33,  also 
indicates  the  quality  of  the  SNR  comparison  Is  good.  The  only  points 
that  lies  outside  the  90  percent  confidence  Interval  are  associated 
with  the  two  point  scatterer  returns. 

Figure  35  shows  the  SNR  comparison  for  a  normalized  filter 
bandwidth  of  0.25.  The  comparison  of  simulation  and  prediction  prior 
to  the  first  point  scatterer  return  shows  that  the  correlated 
scattering  levels  are  3-4  dB  greater  than  the  predicted  uncorrelated 
values,  with  the  greatest  difference  4.3  dB,  at  T/T$  =>  1.5.  Beyond 
the  first  point  scatterer  return,  the  comparison  between  simulated  and 
predicted  levels  appears  to  be  reasonably  good.  The  added  energy  due 
to  the  correlated  scattering  levels  is  still  10  dB  below  the  SNR 
produced  with  the  point  scatterers  returns,  and  therefore  does  not 
appear  to  significantly  add  to  the  overall  SNR. 

The  effects  of  the  correlated  scatter  for  this  filter  bandwidth 
also  appear  in  Figure  36,  normalized  SNR  versus  T/T$.  The  data 
prior  to  the  point  reflection  return  at  T/T$  «=  3.0  is  roughly  3  dB 
above  the  predicted  levels,  and  well  outside  the  90  percent  confidence 
Interval.  The  data  after  the  first  point  scatterer  return  displays 
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more  point  to  point  variation  than  previous  data  sets,  but  no 
systematic  error  Is  apparent.  If  the  correlated  scatterer  affected 
this  portion  of  the  data,  a  similar,  but  smaller  positive  offset  of 
the  data  might  be  expected.  Since  this  Is  not  evident.  It  appears  the 
gain  due  to  correlated  scattering  Is  still  Insignificant  compared  to 
the  point  reflector  returns. 

For  a  filter  bandwidth  of  0.1,  Figure  37,  the  simulated 
correlated  scattering  values  prior  to  the  first  highlight  return  are 
roughly  7  dB  higher  than  predicted  levels.  Beyond  the  start  of  the 
first  point  scatterer  return,  the  SNR  comparison  Is  still  good,  but 
the  simulated  data  shows  a  slight,  ~0.5  dB,  positive  bias. 

The  comparison  for  a  narrow  filter  bandwidth,  0.01,  Is  shown  In 
Figure  38.  The  simulated  level  prior  to  the  first  point  scatterer 
reflection  are  greater  than  the  predicted  levels  Including  the  point 
scatterer  returns.  As  a  result,  the  correlated  continuous  scattering 
Is  the  dominate  scattering  process  at  all  values  of  T/Ts.  Note  the 
simulated  SNR  curve  rises  only  4-5  dB  after  the  point  scattering 
returns,  an  indication  that  the  correlated  scattering  and  highlight 
scattering  energies  are  roughly  equivalent.  The  SNR  gain  prior  to  the 
point  scatterers  Is  roughly  17  dB,  as  listed  In  Table  3. 


CHAPTER  6 


CONCLUSIONS 


6.1  Summary  and  Conclusions 

6.1.1  Models.  In  this  thesis,  results  for  an  analytic  model  and 
computer  simulation  have  been  compared  under  a  variety  of  conditions. 
The  model  assumes  that  a  transmit  pulse  of  duration  T$  propagates 
through  a  lossless  Isotropic  medium  and  Is  reflected  from  a  randomly 
rough,  range  spread  scatterer.  The  scatterer  Is  modeled  as  the  sum  of 
two  components.  The  first  Is  a  continuous  line  scatterer,  and  the 
second,  two  large  amplitude  point  reflectors,  whose  locations  are 
fixed  along  the  length  of  the  continuous  scatterer.  The  reflected 
signal  Is  the  convolution  of  the  transmit  pulse  and  a  random  scatterer 
whose  average  reflection  characteristics  are  described  by  the  target 
scattering  function,  plus  white,  6auss1an  noise. 

The  signal  Is  processed  using  a  complex  crosscorrelatlon 
receiver,  and  the  output  Is  examined  as  a  function  of  Increasing 
receiver  Integration  time,  T.  Integration  times  from  T/T$  «  1.0, 
Integration  time  equal  to  a  transmit  pulse  length,  to  T/Ts  *  11.5,  a 
time  that  Is  slightly  greater  than  the  duration  of  the  reflected 
signal  are  examined.  By  assuming  closed-form  expressions  for  the 
scattering  function  and  transmit  and  receiver  processing  waveforms, 
the  output  of  the  receiver  may  be  expressed  In  terms  of  the  scattering 
function  and  the  receiver  crossambiguity  function.  The  mean  receiver 


output  Is  written  In  terns  of  the  variances  of  the  various  target 
components  and  the  ratio  T/Ts,  predicting  the  receiver  output  signal 
energy  to  noise  energy  ratio. 


6.1.2  Simulation.  These  models  are  simulated  on  a  computer  to 
test  the  validity  of  the  SNR  prediction  equation.  The  continuous 
scatterer  1$  simulated  as  a  set  of  equl-spaced,  uncorrelated,  complex, 
Gaussian,  point  reflectors,  each  of  zero  mean  and  equal  variance. 

This  model  Is  taken  from  VanTrees  for  a  range  spread  scattering 
function.  The  point  scatterer  model  Is  the  slowly  fluctuating  point 
target,  also  from  VanTrees.  The  point  reflectors  are  simulated  as  two 
uncorrelated,  complex,  Gaussian  numbers,  whose  variance  Is  large 
compared  to  the  variance  of  the  continuous  scattering  components.  The 
average  output  of  the  correlation  receiver  Is  computed  by  averaging 
100  realizations  of  the  signal  and  passing  It  through  the  receiver, 
and  separately  100  noise  sequences,  for  constant  scattering  strength 
parameters.  These  two  results  are  then  divided  to  obtain  the 
statistic  to  compare  with  the  prediction  equation. 

6.1.3  Results  and  Errors.  Agreement  between  simulation  and 
prediction  Is  quite  good,  with  the  mean  error  less  than  0.5  dB.  The 
quality  of  agreement  Is  Independent  of  the  absolute  SNR,  the  relative 
scattering  strengths  of  the  scattering  components,  and  Is  generally 
Independent  of  the  receiver  Integration  time.  It  Is  apparent  from 
these  comparisons  that  the  prediction  equation  accurately  predicts  the 
receiver  output  for  all  the  scattering  strength  and  T/Ts  ratios 


examined.  This  further  validates  the  analysis  of  the  model  using 
scattering  and  ambiguity  functions,  allowing  this  technique  to  be 
extended  to  more  complex  scatterers  and  waveforms.  However,  two  types 
of  disagreement  were  Identified,  and  can  be  traced  to  either  a 
transition  between  the  dominate  scattering  components  or  to  the 
effects  of  averaging  over  a  finite  set  of  simulation  values. 

The  first  type  of  comparison  difference  occurs  at  receiver 
Integration  times  that  Include  the  onset  of  the  first  highlight 
reflection.  These  times,  at  values  of  T/Tg  ■  3.0  and  3.5  correspond 
to  processing  the  reflected  signal  on  the  start  of  the  first  point 
scatterer  reflection.  For  T/Ts  -  3.0,  only  a  single  sample  of  the 
large  amplitude  reflection  Is  Included  with  the  continuous  scattering 
return.  For  T/Ts  -  3.5,  half  the  first  point  scattering  return  Is 
Included.  The  comparison  difference  for  data  at  T/T$  ■  3.5  Is 
generally  2-3  dB  greater  than  the  overall  mean  comparison  difference. 
Since  a  small  number  of  large  values  are  being  averages,  when  compared 
to  either  greater  or  smaller  values  of  T/T$  are  averaged.  It  Is 
expected  that  the  variance  In  the  estimate  will  be  greater.  This 
behavior  Is  not  seen  In  the  comparison  at  T/T$  «  3.0  because  the 
energy  In  the  single  point  scatterer  sample  Included  by  the  receiver 
Is  small  compared  to  the  energy  In  the  300  continuous  scattering 
samples,  and  does  not  effect  the  estimate  of  signal  energy,  as  does  50 
samples  for  T/T  -  3.5. 


An  Independent  confirmation  of  this  conclusion  comes  by 
Increasing  the  number  of  model  realizations  averaged  to  obtain  the 
simulation  comparison  value.  While  this  Is  cost  prohibitive  for  all 
combinations  of  scatterers  and  receiver  Integration  times,  comparisons 
of  1000  and  10000  averages  show  the  comparison  error  at  T/T$  -  3.5 
declines  dramatically.  Indeed,  for  large  numbers  of  signals  In  the 
average,  all  the  simulation  data  approaches  the  predicted  values.  It 
Is  concluded  that  the  differences  seen  between  simulation  and 
prediction  outside  the  region  near  T/Ts  ■  3.0,  the  second  type  of 
error  observed,  can  be  reduced  by  Increasing  the  number  of  model 
realizations  used  to  compute  the  average  signal  and  noise  energies. 

6.1.4  Correlated  Scattering.  A  modification  of  the  continuous 
scattering  model  was  made  by  filtering  the  continuous  scattering 
sequence  to  produce  a  sequence  with  correlation  between  the  continuous 
scattering  components.  The  scattered  signal  Is  processed  as  described 
above,  and  the  averaged  result  compared  to  the  prediction  for 
uncorrelated  scattering. 

For  the  correlated  continuous  scatterer  only,  the  quality  of  the 
comparison  degrades  with  decreasing  filter  bandwidth.  In  all  cases, 
the  simulated  SNR  Increased  above  the  prediction  with  decreasing 
bandwidth.  The  largest  mean  difference  observed  Is  17  dB,  for  a 
normalized  filter  bandwidth  of  0.01.  This  Increase  In  energy  over 


the  uncorrelated  value  Is  termed  coherent  SNR  gain.  It  appears  the 
correlation  between  the  continuous  scattering  components  adds  coherent 
energy  to  the  reflected  signal,  more  like  a  deterministic  spread 
component. 


This  effect  Is  also  observed  with  the  composite  scatt<~  r 
comparison.  However,  the  large  point  scatterer  returns  can  swamp  out 
the  coherent  energy,  except  for  very  small  filter  bandwldths.  For 
mildly  correlated  scatterer  the  prediction  and  simulation  disagree 
prior  to  the  first  highlight  return,  but  the  energy  In  the  point 
scatterer  returns  Is  large,  and  beyond  T/Ts  ■  3.5,  the  simulation 
and  prediction  agree  quite  well.  With  very  correlated  scattering,  the 
coherent  energy  dominants,  and  the  point  scatterer  returns  are  barely 
noticeable  for  all  values  of  T/Ts.  In  general,  a  filter  bandwidth 
of  0.05  or  less  was  necessary  for  the  coherent  scatter  to  dominate, 
using  a  point  scattering  to  continuous  scattering  strength  ratio  of 


The  encouraging  results  of  this  thesis  Imply  the  ambiguity 
function,  scattering  function  approach  can  be  applied  to  a  wide 
variety  of  scatterers  and  receiver  waveforms.  This  Is  not  a  new 
conclusion  or  result,  rather  a  further  demonstration  of  fact.  The 
simulation  results  and  predictions  used  In  this  thesis  In  particular 
suggest  further  efforts  In  two  areas. 


The  first  area  Is  to  Increase  the  complexity  of  the  scattering 
model.  Straightforward  modifications  might  Include  a  time  varying 
continuous  scattering  amplitude,  several  point  scatters  with  random 
locations.  The  prediction  equation  can  be  similarly  modified,  subject 
to  the  restriction  that  the  modification  can  be  expressed  analytically. 
In  either  a  deterministic  or  statistical  sense.  Adding  Doppler 
spreading  Is  also  possible  In  the  model,  but  more  difficult  to  Include 
In  the  analytic  prediction. 

The  other  set  of  potential  modifications  to  the  model  Included  an 
Investigation  of  the  effects  of  various  processing  waveforms.  This  Is 
a  simple  change  In  the  model,  but  requires  a  reworking  of  the 
prediction  equation  to  Include  the  new  crossambiguity  function.  The 
value  of  this  type  of  modification  allows  evaluation  of  very  complex 
waveforms  with  a  simple,  well-behaved,  scatterer. 
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APPENDIX  A 


DERIVATION  OF  THE  CROSSAMBIGUITY  FUNCTION 

The  general  ambiguity  function  Is  defined  as 


A(  t  ,  «')  - 


where 
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U  *  «  -  «Q  , 

where  rQ  and  «q  are  the  delay  and  doppler  of  the  actual  scattering 
function  and  t  and  w  are  the  delay  and  Doppler  assumed  by  the  receiver. 


ACVtu1)  Is  a  measure  of  the  mean  square  deviation  of  the 
transmit  waveform  with  a  delayed  and  shifted  version  of  Itself  [22]. 
The  function  A(t' ,»')  Is  referred  to  as  the  autoambiguity  function  of 
the  signal  f(t).  This  result  may  be  extended  to  two  different 
waveforms.  The  crossambiguity  function  of  f(t)  and  g(t)  Is  then 
defined  as. 
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and  Is  a  measure  of  the  mean  square  error  of  one  signal.  In  our  case 
the  transmit  signal,  with  a  delayed  and  shifted  version  a  second 
signal,  referred  to  as  the  receiver  processing  signal. 

T(t)  and  ^(t)  are  assumed  to  be  of  the  form. 


where  T  >  T  . 
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Equations  3  and  4  Imply  the  envelope  of  the  transmit  signal  Is 
assumed  to  fall  completely  within  the  duration  of  the  processing 
signal  and  T$  <  T.  Substitution  of  these  signals  Into  Equation  A-2 
yields 
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APPENDIX  B 

GENERATION  OF  RANDOM  NUMBERS 


B.l  Generation  of  Uniform  Random  Numbers 

The  target  scattering  simulation  uses  the  operating  system 
supplied  uniform  random  number  generator,  RAN,  to  compute  uniform 
numbers  In  the  range  [0,1]  [23].  The  first  call  to  the  generator 
requires  a  unique  Integer  seed,  which  should  also  be  large  and  odd. 
This  seed  Initializes  the  generator,  and  Is  calculated  by  doubling  the 
value  of  the  computer  system  time  clock.  In  seconds,  and  adding  one. 
Twenty  additional  calls  are  then  made  to  RAN,  to  produce  the  first 
seed  value  used  for  the  scattering  simulation.  Since  the  simulation 
data  was  computed  over  a  number  of  days,  none  of  the  original  seeds 
were  Identical,  and  It  Is  reasonable  to  assume  the  sequences  of  random 
numbers  generated  used  In  the  simulation  are  statistically  Independent 
from  each  other. 

B.2  Generation  of  Gaussian  Numbers 

Each  Individual  scatterer  Is  modeled  as  a  sample  function  from  a 
complex  6au$s1an  process.  To  generate  these  scatterers  In  simulation, 
two  real.  Independent  Gaussian  numbers  are  required.  Each  pair  of 
numbers  Is  obtained  by  successive  calls  to  the  Gaussian  number 
Generator  GGNHL,  a  subroutine  In  the  IMSL  Mathematical  Library  [24]. 
The  generator  output  was  tested  for  two  statistical  properties; 


distribution  of  generator  outputs,  and  correlation  between  sequences 
generated  with  different  Initialization  seeds. 

Figures  39  and  40  are  histogram  plots  for  G6NML  outputs  for  100 
and  1000  sample  sets,  respectively.  While  the  theoretical  values  of 
the  mean  and  standard  deviation  are  0  and  1,  the  histogram  of  the  100 
sample  set.  Figure  3  has  a  mean  of  0.04  and  a  standard  deviation  of 
0.92.  A  chi-squared  test  was  performed  on  this  set,  and  the 
assumption  of  6auss1an  distribution  was  satisfied  with  80  percent 
confidence  [25].  For  the  1000  sample  set,  the  mean  and  standard 
deviation  are  0.01  and  1.01,  very  close  to  the  values  of  the 
continuous  distribution.  The  confidence  that  this  sequence  Is 
6auss1an  distributed  Is  greater  than  95  percent.  For  sets  of  greater 
than  10000  samples,  the  mean  and  standard  deviation  are  equal  to  the 
continuous  distribution  values,  and  on  the  basis  of  these  tests.  It  Is 
concluded  that  the  random  number  generator  GGNML  does  produce 
sequences  of  numbers  that  are  Gaussian  distributed. 

The  second  test  applied  to  these  sequences  gives  a  measure  of  the 
Independence  between  samples  generated  with  different  Initial  seed 
values.  Figure  41  Is  a  representative  plot  of  the  autocorrelation  of 
a  1024  sample  set.  The  ordinate  of  the  plot  Is  correlation  magnitude, 
normalized  between  -1  and  1.  The  absclcca  Is  delay,  In  units  of 
samples.  The  correlation  magnitude  Is  less  than  0.1,  beyond  10 
samples  of  delay.  This  Indicates  that  the  sequence  Is  poorly 


SAMPLES 


Autocorrelation  of  random  number  generator 
output  for  a  sample  set  size  of  1024. 


correlated  with  Itself,  a  Indication  that  successive  generator  outputs 
are  uncorrelated. 

A  representative  crosscorrelatlon  between  two  sequences  generated 
with  different  Initial  seed  values  Is  shown  In  Figure  42.  Unlike  the 
autocorrelation  example,  where  the  correlation  magnitude  at  a  delay  of 
zero  samples  Is  unity,  the  crosscorrelatlon  magnitude  Is  low  across 
all  delay.  These  results  Indicate  that  the  Gaussian  number  generator 
outputs  are  reasonably  uncorrelated  with  each  other.  This  result  Is 
expected,  since  the  theoretical  distribution  values  are  defined  to  be 
statistically  Independent  [26].  While  statistical  Independence  of  the 
generator  output  cannot  be  Inferred  on  the  basis  of  correlation,  that 
the  outputs  show  poor  correlation,  both  auto  and  cross.  Is  sufficient 
for  this  thesis  to  validate  the  random  number  generator  output  [27]. 


APPENDIX  C 


CALCULATION  OF  VARIANCE  OF  SIGNAL-TO-NOISE  RATIO 


The  slgnal-to-nolse  ratio  Is  computed  by  evaluating 
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where  and  are  zero  mean  Gaussian  variables.  Therefore  Sk  is 

R  I  2 

with  2  degrees  of  freedom.  Similarly,  Nfc  is  X  with  2  degrees  of 

freedom.  The  sum  of  N  chi-squared  variables  is  also  chi-squared,  with 

2N  degrees  of  freedom. 
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From  Nuttall  [28],  the  PDF  of  gamma  variates  Is 
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and  the  PDF  of  the  ratio  z  -  x/y  for  x  and  y  Independent  Is 
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For  the  case  of  Interest  y  ■  »  •  H,  an  Integer.  Equation  6  may 
be  reduced  to. 
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From  6radshteyn  and  Ryzhlk  [29], 
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and  the  variance. 
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Equation  C-17  Is  Identical  In  form  to  the  variance  of  the  F 
distribution  with  2N  degrees  of  freedom  [31].  It  Is  therefore 
concluded  that  variance  of  the  simulated  slgnal-to-nolse,  normalized 
by  the  ratio  of  the  noise  variance  to  signal  variance  Is  F 
distributed.  Confidence  limits  may  now  be  computed  by  calculating 
F(2N,2N)  for  a  specific  a. 
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